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D scriplion 

Field of Invention 

5 This invention relates to directional coupler filters and more particularly to such filters wherein the spacing between 

the waveguides is tailored to maximize sidelobe suppression and to minimize bandwidth. 

Background of Invention 

10 The potential communication capacity of optical fibers operating in the low loss wavelength windows of 1 .3 um 

and 1.5 jim is in the order of tens of Terahertz. The practical utilization of this bandwidth may be realized through the 
use of wavelength division multiplexing (WDM). In this scheme the spectral range is subdivided and allocated to different 
carrier wavelengths (channels) which are multiplexed onto the same fiber. The frequency bandwidth that an individual 
channel occupies depends on a number of factors, including the impressed modulation bandwidth, margins toaccom- 

15 modate for carrier frequency drift and carrier frequency uncertainty, and also to reduce possible cross-talk between 
channels. 

Although the isolated optical fiber may inherently have tremendous information carrying capacity, the overall optical 
communication link may be significantly restricted in bandwidth. These restrictions may result from the limited optical 
amplifier spectral windows, the availability of lasing sources and their tuning ranges, and also filter tuning ranges. 
20 Hence, to achieve efficient use of bandwidth requires that the available communications windows be densely filled with 
multiplexed channels. At the output of such a system, filters are needed to separate the wavelength channels. The 
performance of these wavelength filters in their ability to reject out of band signals, critically determines channel spacing 
and hence channel density. 

The type of filters of interest here operate through a wavelength dependent exchange of power between two 

25 waveguide modes. It is well known that two waveguides placed in close proximity may exchange power through their 
, evanescent fields, which penetrate the guiding layer of the other waveguide. This power exchange occurs continuously 
along the propagation direction, with a rate determined by the inter-waveguide spacing and the degree of velocity- 
matching of the two modes. If the velocities of the two waveguide modes are identical, the situation is termed 'syn- 
chronous' or 'phase matched', and the power coupled into one of the waveguides accumulates constructively. Complete 

30 power exchange is then possible and occurs at a characteristic coupling length L c which is determined by the structure 
of the device. If the modes propagate at different velocities, then this condition is termed 'non-phase matched' or 
'asynchronous'. The power in the coupled waveguide accumulates with a phase error, leading to incomplete power 
transfer in this case. The larger the phase-mismatch, the faster the phase-error accumulates, which results in less 
power being transferred to the coupled waveguide. 

35 FIGURE 1 a shows two typical waveguides placed in parallel, a configuration known as the directional coupler. The 

input power is initially launched into waveguide 1 (the 'input' guide), and the output is extracted from guide 2 (the 
'coupled' guide). FIGURE 1b shows the power in waveguide 2 (the 'coupled' guide) as a function of device length for 
two cases of phase-matching. The solid curve represents a synchronous case, with complete power exchange occur- 
ring at L c = 5mm. The dashed curve represents coupling between asynchronous modes. 

40 Wavelength selectivity in the directional coupler occurs through differential velocity dispersion. At the design wave- 

length, the velocities of the two modes are equal. As the wavelength is changed or 'detuned', the mode velocities 
necessarily change. For filtering action however, it is critical that the difference between these mode velocities changes, 
i.e., a differential velocity dispersion is required. The rate of change of differential velocity with respect to wavelength 
is the primary factor in determining filter bandwidth. This rate is a function of material type and waveguide structure. 

45 FIGURE 2 shows the filter response of a parallel directional coupler. The abscissa is in terms of a normalized detuning 
factor, Ap, which is a measure of the velocity difference between the coupled modes. This axis can be converted into 
an actual wavelength scale when the relationship between A(3 and wavelength \, is established for a particular device. 
The ordinate is the power in the coupled waveguide, in logarithmic scale, for a device of fixed length. The normalized 
half-power bandwidth is 8.4 radians (rad), and the maximum sidelobe level is at -9.3 db. 

so For optical communications purposes, a sidelobe level of -9.3 db is too large, since it would represent a significant 

cross-talk to an adjacent wavelength channel, if these channels were spaced by the width of the passband. If it is 
required that the cross talk in an adjacent channel be less than -9.3 db, the spacing between adjacent channels in the 
wavelength domain must be made much larger than the main passband width. Since the sidelobe levels decrease at 
a slow rate with detuning (see FIGURE 2), the channels must be widely separated. Hence, a severe penalty is paid in 

55 , terms of channel density, and hence information carrying capacity, for the price of low cross-talk. It is very desirable 
then to identify some degre of freedom which may be used to improve filter response. 

The degre of freedom most commonly used in directional couplers is a modulation of the interaction strength of 
the two coupled modes. This may be achieved for example, by modulating the inter-waveguide separation in the di- 
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rectional coupler. In many other branches of optics and physics this process is known as 'apodizationV In waveguide 
theory it is referred to as tapering'. 

Directional coupler devices are commonly modeled through a set of coupled differential equations such as: 

5 ^ = .j kG ^ z A 2 Equation (1a) 

^ = -jke!^ z A^ Equation (1b) 

10 In (1a) and (1b) Ay and A% represent the amplitudes in waveguides 1 and 2. AP is the detuning constant and k is the 
interaction strength, k depends on the waveguide structure and is strongly influenced by the separation between ad- 
jacent waveguides. The physical origin of this coupling may be due to the interaction of evanescent fields in a uniform 
coupler, (as in FIGURE 1a), or coherent scattering in a grating-assisted coupler. 

By varying the interaction strength k along the directional coupler, the spectral response of the device may be 

is improved. The physical origin of this improvement is in the interferometric nature of the coupling process: at every 
position along the coupler, power is being transferred from the input waveguide to the coupled waveguide. The total 
power in the coupled waveguide at some point then, is an interferometric sum of all the power coupled into the 
waveguide prior to that point. That is, a sum including relative phase delays. By adjusting the interaction strength k 
along the waveguide, one dictates the rate of power transfer at each position, abng with its phase relationship to the 

20 total coupled power. By judicious choice of the coupling taper shape A(z), it is theoretically possible to generate any 
(passive) response. 

Howtocalculate the taper shape k(z) for a desired response has been a longstanding unanswered design question. 

The original proposal, suggested in 1978 by Alferness et al, IEEE J. Quantum Physics, QE-14. No. 11, 1978, pp. 

843-847, was based on an approximate Fourier transform relation which gave a few promising shapes. This proved 
25 to be a useful guideline in improving actual device response (Alferness, Applied Physics Letters, Vol. 35, No. 3, 1979, 

pp. 260-262), but thus far has been unsuccessful in yielding the sidelobe suppression required in communication 

systems. To date no improvements have been advanced, and the same shapes suggested in 1978 continue to be the 

only ones analyzed (see, for example, H. Sakata, Optical Letters, Vol. 17, No. 7, 1992, pp. 463-465). 

The goal of filter design is to solve for the interaction function k(z) of (1 ), given some desired output response as 
30 a function of wavelength 4g(X). However, when k(z) is non-constant, the set of coupled equations in (1 ) has no analytic 

solution in general. Hence, filter design is currently guided by a set of approximate solutions. The most important of 

these approximate solutions is the Fourier transform relation, given by: 

A>(AP) = J k(z) e^ M tfz Equation (2) 

35 in Equation (2), /^(Afi) is the amplitude in the output or coupled waveguide as a function of detuning AP, (which may 
be related to the actual wavelength X). Because Equation (2) represents a Fourier transform relation between k(z) in 
the spatial domain and /^(AP) in the wavelength domain, the principle of duality may be used. That is, given a desired 
4 ? (AP), k(z) is found by the inverse Fourier transform. This approximation is valid for small coupling values, and does 
not extend to describe the critical region of the main passband and first few sidelobes. No analytic solution currently 

40 treats the important region around the central wavelength. 

Summary of the Invention 

The present invention seeks to produce a directional coupler filter having a specific sidelobe level. 
45 The present invention also seeks to produce a directional coupler filter having a specific sidelobe level and a 

minimum spectral bandwidth for the specific sidelobe level. 

The present invention seeks to produce the aforementioned directional coupler filter by specifying an interaction 
strength function between waveguides. 

The present invention seeks to produce the specified interaction strength by defining the shape function between 
50 the waveguides. 

Therefore in accordance with a first aspect of the present invention there is provided a directional coupler filter 
comprising a first waveguide for receiving an optical signal having at least one wavelength; a second waveguide in a 
spaced relationship to the first waveguide so that the wavelength may be selectively coupled from the first waveguide 
to the second. The spaced relationships is tailored so as to suppress the spectral sidelobes of the coupled wavelength 
55 to a level down to -75 db. 

In a preferred embodiment the spectral bandwidth is a minimum for the selected sidelobe.. 

In accordance with a second aspect of the present invention there is provided a method of suppressing sidelobe 
levels in an optical directional coupler f ilt r, th filter having a first waveguide for receiving an optical signal with at least 
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one wavelength and a second waveguide in spaced relationship to the first waveguide, the wavelength being selectively 
coupled to the second waveguide, the method comprising tailoring the interactive strength of th first and second 
waveguides according to formula: 

s k(z) = L 0 (z) + SL, (z) + S%(z) + S 3 L 3 (z) 

wherein 

k{z) - interactive strength, 

S= desired sidelobe level, in absolute 6B, (IdSI), and 

10 Lj = LpLpL&Lg,... are functions of propagation distance (z). 

Brief Description of Drawings 

The invention will now be described in greater detail with reference to the appended drawings wherein: 

15 

FIGURE 1a is a plan view of a directional coupler having parallel (non-tapered) waveguides; 

FIGURE 1 b shows the power in the coupled waveguide as a function of propagation distance for synchronous and 

asynchronous designs; 

FIGURE 2 illustrates graphically the spectral response of a parallel coupler; 
20 FIGURE 3 shows the filter bandwidth as a function of the highest sidelobe level; 

FIGURE 4a shows coupling strength as a function of normalized propagation distance for an optimum coupler with 
-60 db sidelobes; 

FIGURE 4b is the spectral response for the filter of FIGURE 4a; 

FIGURE 4c is the spectral response of a physical device simulated with the beam propagation numerical method 
25 (BPM) and using the shape of FIGURE 4a; 

FIGURE 4d illustrates physical separation between two adjacent waveguides for the BPM simulation of FIGURE 4c; 
FIGURE 5 shows taper shapes of -40 db, -50 db and -75 db couplers generated in accordance with the formula 
of the present invention; 

FIGURES 6a, 6b and 6c show spectral responses of the couplers of FIGURE 5 - for -40 db, -50 db and -75 db 
30 respectively; 

FIGURE 7 illustrates the bandwidth/sidelobe level relationship for a coupler designed in accordance with the 
present invention and the theoretical optimum; 

FIGURE 8 is a plan view of a practical embodiment of the present invention; 
FIGURE 9 is a plan view of the coupler of FIGURE 8 having gratings to assist the coupling; and 
35 FIGURE 10 is a plan view of the device of FIGURE 8 having electrodes for electrooptic tuning. 

Detailed Description of the Invention 

As stated the objective herein is to produce filters with a specific sidelobe level. The response of the coupler in 
40. between sidelobes is unimportant, since by definition the sidelobes represent local extrema. If we have some coupler 
shape f^fz), then the following error vector is defined: 

£(*„) = [£c{k 0 ),eM- £ M]- Equation (3) 

45 The components z- p {M ...m] describe the error between the desired sidelobe level and the current sidelobe level 

for shape Aq, for each of the m sidelobes considered: e,- = (a ; - p/), where a,- is the desired sidelobe level and p, is the 
current sidelobe level for sidelobe /. 

We let the arbitrary coupling coefficient be represented as A(z) = Z N B n f n (z) t where the /„ is an arbitrary set of 
functions. The unknowns are the coefficients a n which are written in a vector, K = (a 1( %,...a w ). Given an initial state 

50 £ 0 , an improvement to K 0 denoted by £* = K 0 + A £ is found by evaluating the incremental improvement vector Ate, 

Ak = oJ T (jJ T y l E{k 0 ) Equation (4) 

55 where J Fj =^ is the Jacobian, and a is a scaling parameter used to decelerate the convergence. Equation (4) is 

repeatedly applied until a satisfactory optimization is achieved. Hence the optimum value is analytically given by the 
function k(z) s 2 a n f n (z). 
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Achieving th desired sidelobe level is the first criteria in filter synth sis. The second criteria is to obtain the nar- 
rowest bandwidth while still maintaining the maximum tolerable sidelob level. This second criteria is met when all 
sidelobe levels are at the maximum tolerable level. This may be proved rigorously in the following way. Define the 
coupler bandwidth B, as a function of the levels of all the sidelobes: B s BfS^S^S^...) where S y - is the absolute value 
of sidelobe / Evaluate the gradient of B with respect to the vector^ = (s 1 ,s 2l s 3l ...) > 



Equation (5) 



\7v*(s) is numerically shown to be positive definite, indicating that decreasing any sidelobe has the effect of increasing 
the bandwidth B. Hence for a desired maximum sidelobe level x, all sidelobes must be below x by definition, but to 
achieve the minimum bandwidth, they should only be infinitesimally below x. 

Using the above method, a coupler with arbitrarily low sidelobes may be theoretically generated. Before actual 
shapes are presented for various couplers, the overall achievements are compared to currently known taper functions. 
In FIGURE 3, the trade off between the desired sidelobe level (abscissa) and the bandwidth at that level (ordinate) is 
shown. The solid curve represents the theoretical results derived here. This represents the narrowest possible band- 
width at a desired sidelobe level, a coupler cannot have a response that lies below this line. For comparison, the 
response of a few of the well known shapes are plotted as labeled points. The best known couplers have a maximum 
sidelobe level in the range of -40 db to -45 db. On the other hand, the method presented here may generate couplers 
with sidelobe levels several orders of magnitude lower. 

A representative result of using Equations (3) and (4) for filter synthesis is given in FIGURE 4a, which shows the 
optimum shape of k(z) for a coupler designed for -60 db sidelobes. The abscissa is in normalized length units. The 
actual physical length of the device is obtained by multiplying the abscissa in FIGURE 4a by the coupling length 
of the specific device. The ordinate is the strength of k(z) in normalized units. k(z) may be translated into either inter- 
waveguide spacing or grating strength, in a specific uniform coupler or grating-assisted coupler, respectively. FIGURE 
4b shows the theoretical spectral response of the taper shape of FIGURE 4a. All the sidelobes are at -60 db, in order 
to achieve the minimum bandwidth for this maximum tolerable sidelobe level. FIGURE 4c shows the response of a 
numerically simulated physical device. The simulation was performed using the Beam Propagation Method (BPM). 
This figure shows the very close correspondence that may be achieved in practice. The actual shape of the coupler is 
shown in FIGURE 4d, which depicts the total separation between two adjacent waveguides (in urn) as a function of 
propagation distance (in cm). 

A practically useful range for filter sidelobe design is the range of -40 db to -75 db. An approximation to the optimum 
shapes which cover this range is given by the design formula 

k(z) = L^z^fz^L^^fz^L^z), 
where S is the desired sidelobe level in absolute decibels and the functions Lfz) are given by 



Equation (6) 



40 



W = £^ y cos((2y-l)TO) T (-0.5<z<0.5) 



Equation (7) 



45 



50 



55 



The constant coefficients by are given in TABLE 1 . 

TABLE 1 





i | 


b u 


1 


2 


3 


4 


5 


6 


b 0j 


1.96526 


-0.175916 


1.03664 


-0.514003 


0.468375 


-0.335236 


b U 


2.65083e-2 


-7.15423e-3 


-6.89689e-2 


3.01025e-2 


-2.79768e-2 


1.99425e-2 


b 2J 


-2.82597e-4 


1.24012e-3 


1.68102e-3 


-6.50960e-4 


6.26496e-4 


-4.46475e-4 


b 3J 


8.61456e-7 


-1.92417e-5 


-1.73664e-5 


6.26998e-6 


-6.23891 e-6 


4.45589e-6 


b 4J 


2.901 13e-9 


9.51447e-8 


6.91843e-8 


-2.21879 -8 


2.34583e-8 


-1.66805e-8 



5 



EP 0 726 476 A1 



The normalized propagation distance is z, which ranges over {-0.5<z<0.5}. For a real devic , the physical length scales 
as ZszLp where Zis the physical length, z is the normalized length and L c is the length of the interaction r gion in the 
device. By substituting for a value of S in the range 40 to 75, Equation (6) gives the required interaction strength. The 
shapes for -40 do, -50 db, and -75 db sidelobe suppressed couplers are shown in FIGURE 5. The corresponding 
s spectral responses are shown in FIGURES 6a-6c. The trade-off in bandwidth versus sidelobe level for couplers de- 
signed by Equation (6) is shown in FIGURE 7. The solid curve represents the results of Equation (6) while the dashed 
curve represents the ideal cases (as shown in FIGURE 3). 

FIGURE 8 illustrates, generally, one practical embodiment of the present invention. It is to be understood that the 
spacing between waveguides must be such that the interactive strength satisfies Equation (6). This may be achieved 
10 by curving both waveguides as shown in FIGURE 8 or by curving only one while the second remains straight. 

It is known that the interaction strength between the two waveguides depends exponentially on their separation, 

K = K Q exp{-yd) Equation (8) 

where dis the separation between the two waveguides (edge to edge). In Equation (8), Kq and yare constant coefficients 
15 that need to be determined for a specific device. By selecting two values of separation d, and d 2 , two corresponding 
values of coupling strength, /C, and K 2 may be determined by well known methods (see Haus et al, IEEE J. Lightwave 
Technology, Vol. LT-5, No. 1 , pp. 16-23, 1 987). Hence we may solve for the coefficients of Kq and y. 
Equation (8) may be rearranged to write d in terms of K: 

K L 

20 d= \ ln{ -h^ ) Equation (9) 

K L C is the normalized coupling strength of Equation (6), (that is, k(z) = K L C ). Hence the physical separation d, 
given in terms of the normalized coupling strength k(z) is 

K L 

25 d(Z) = * ln( ) Equation (10) 

where Z is the physical distance, L c is the coupler's interaction length, z is the normalized distance given by z = 
and k(z) = k(Z/L^ is the function given in Equation (6). 

These values can then be used to plot dvsZox the relationship can be incorporated into mask layout design 
30 software to generate the necessary pattern to produce the waveguides. 

The practical embodiment of the invention can be implemented in semiconductor material such as 111 - V or ll-VI 
alloy compounds as well as a Si/Ge system. The configuration is also applicable in silica, glass, polymers and pho- 
torefractive materials such as lithium niobate. 

While particular examples of the invention have been described it will be apparent to one skilled in the art that 
55 variations and alternatives are possible. Such variations include both lateral and vertical coupling of the waveguides. 
Also included are waveguides in which one or both are provided with gratings 4, as shown in FIGURE 9. It is also 
contemplated that the waveguides can be provided with electrode means 6 (FIGURE 10) in order to effect electro-optic 
tuning. It is to be understood, however, that such variations and alternatives fall within the scope of the invention as 
defined by the appended claims. 

Claims 

1. A directional coupler filter comprising: 

a first waveguide for receiving an optical signal having at least one predetermined wavelength; 
a second waveguide in a spaced relationship to said first waveguide so that said predetermined wavelength 
is selectively coupled from said first waveguide to said second waveguide wherein said spaced relationship 
is tailored so as to suppress spectral sidelobes of said coupled wavelength to a level down to -75 db. 

2. A directional coupler filter as defined in claim 1 , said spaced relationship determining the interactive strength be- 
tween the first and second waveguides. 

3. A directional coupler filter as defined in claim 2, the interactive strength being related to sidelobe level and prop- 
55 agation distance. 

4. A directional coupler filter as defined in claim 3, fabricated in accordance with the formula 



45 



so 
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k(z) = L 0 (z)+SL, (z)^L 2 (z)^L z (z) 

where 

5 k(z)= interactive strength, 

S sidelobe lev I in IdBI, and . 

L^L^Lj.Lj,... are functions of the propagation distance (z). 

5. A directional coupler filter as defined in claim 4, fabricated in a semiconductor material. 

10 

6. A directional coupler filter as defined in claim 5, said semiconductor material being a lll-V alloy compound. 

7. A directional coupler filter as defined in claim 5, said semiconductor material being a ll-VI alloy compound. 
15 8. A directional coupler filter as defined in claim 4, fabricated in a glass material. 

9. A directional coupler filter as defined in claim 4, fabricated in a polymer material. 

10. A directional coupler filter as defined in claim 4, fabricated in a photorefractive material. 

20 

11. A directional coupler fitter as defined in claim 10, said photorefractive material being lithium niobate. 

12. A directional coupler filter as defined in claim 1, at least one of said waveguides having gratings to assist the 
coupling. 

25 

13. A directional coupler filter as defined in claim 1 , at least one of said waveguides having electrode means in order 
to supply injection current thereto. 

14. A directional coupler filter as defined in claim 4, said sidelobe level dependent on the spectral bandwidth of said 
30 coupled wavelength. 

. 15. A directional coupler as defined in claim 14, said spectral bandwidth selected to be a minimum. 

16. A method of suppressing the sidelobe levels in an optical directional coupler filter having a first waveguide to 
35 receive an optical signal having at least one predetermined wavelength and a second waveguide in spaced rela- 

tionship to said first waveguide, said predetermined wavelength being selectively coupled to said second 
waveguide, the method comprising: 

tailoring the spaced relationship between said first and second waveguides so as to provide an interactive 
strength according to the formula: 

m = lofcH-a, (zh^L^z)^L z (z) 

wherein 

k(z)~ normalized interactive strength, 

45 S= desired sidelobe level in IdBI, and 

Lq.L) - are functions of propagation distance (z) 



so 
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